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Abstract 

We study the dynamical exponents dw and ds for a particle diffusing in a 

disordered medium (modeled by a percolation cluster), from the regime of 

extreme disorder (i.e., when the percolation cluster is a fractal at p = pc) to 

the Lorentz gas regime when the cluster has weak disorder at p > pc and the 

leading behavior is standard diffusion. A new technique of relating the velocity 

autocorrelation function and the return to the starting point probability to the 

asymptotic spectral properties of the hopping transition probability matrix of 

the diffusing particle is used, and the latter is numerically analyzed using the 

Arnoldi-Saad algorithm. We also present evidence for a new scaling relation 

for the second largest eigenvalue in terms of the size of the cluster, | In AmaxI ~ 

g-dyj/df ^ which provides a very efficient and accurate method of extracting 

the spectral dimension dg where ds = 2df /dw 
05.40.+j, 05.50.+q, 64.60.Fr 
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I. INTRODUCTION 



Random walks have held the interest of physicists, engineers and mathematicians for a 
long time now, partly due to their ability to model a wide variety of problems which can 
be described at some level by transition in states which are governed by probabilistic laws. 
For example, random walks have been used extensively to model the motion of a particle in 
diffusive processes through solids (ordered or disordered), liquid or gas, including transport 
problems like electrical conduction ||T],0 , and in describing certain conformational properties 
of macromolecules [^,0. 

In this paper we present an approach to studying random walks which is an alternative 
to the more common one where the individual random walk is pursued either by simulation 
or in some other way. This is the method of spectral analysis of the corresponding transition 
probability matrix. This method allows us to relate the various dynamical quantities of the 
walk, such as the velocity autocorrelation function, {'v{t) ■ v(0)), the return to the starting 
point probability after t steps, P(t), the mean square displacement, and the acceleration 
autocorrelation function, to the asymptotic spectral properties of the transition probability 
matrix W (defined below). With the development of the various numerical methods of 
diagonalizing large matrices, this approach can serve as a very powerful tool to solve for 
these dynamical quantities. 

In this work, we model a random medium by a percolation cluster, where sites on a 
lattice are occupied with probability p (otherwise empty) and the occupied nearest neighbor 
sites are considered to be connected It is well known that there is a critical threshold 
Pc where the incipient infinite cluster is a disordered fractal . Such a cluster will serve as 
an example of a random medium with a critically correlated disorder. We will also confine 
our attention for the moment to the random walk which takes place on the (critical or 
noncritical) percolation cluster according to the so-called blind or myopic ant rule. The 
blind ant moves to an occupied nearest neighbor in one time step with equal probability of 
1/z (where z is the coordination number of the underlying lattice) and it remains at the 
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same site with the probabihty equal to {z — Zi)/z (where Zi is the number of occupied nearest 
neighbors of that particular site). On the other hand, the myopic ant never stays at the 
same site but always moves to an occupied nearest neighbor with equal probability l/zi at 
each time step 

The transition probability matrix W is defined as the matrix with its ij-ih element equal 
to the probability that a particle at site j hops to site i in the next time step. Both the 
blind and myopic ants can be described by a matrix W which is a Markov matrix P| as each 
element is equal to or greater than zero and the sum of all the elements of each column is 
equal to one. Such a matrix describes a Markov chain, which is a one-dimensional, discrete 
Markov process. 

In this paper we will mainly concentrate on the two important quantities for the random 
walker, (v(t) ■ v(0)) and P{t). At p = Pc blind and myopic ants are known to exhibit a 
nondiffusive leading behavior called anomalous diffusion 0. In this case, the walk is often 
characterized by a walk dimension where its mean square displacement after t steps 
{R{tY) is given by 



The deviation of d^ from 2 is a measure of the nonstandard nature of the random walk. 

In a stationary ensemble (that is, where the correlation functions depend only on the 
relative time difference and not on the absolute times), the second derivative of (R(i:)^) is 
proportional to (v(i(:) • v(0)). So under these conditions (v(t) • v(0)) is normally expected to 
obey the following power law in the long time limit: 



Such a behavior is common for the blind ant and also for the myopic ant on so-called 
nonbipartite lattices at least asymptotically (although the leading behavior in the velocity 
autocorrelation for the myopic ant on bipartite lattices do not derive from that of {R{tY), 
thus providing exceptions to this behavior). A lattice (or a cluster) is called bipartite if it 



(1) 



(v(t) ■ v(o)) ~ 



(2) 
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can be partitioned into two equivalent classes of sites where all connectivity is between sites 
of different classes. 

The long-time behavior of the probability P{t) that a random walker returns to the 
starting point at time t is also generally a power law, where 

p{t) ~ r-^^/^ (3) 

and ds is the spectral dimension. For a standard diffusive behavior, ds is equal to the spatial 
dimensionality d, while on a fractal, it is generally a fractional number different from d. 
Thus one can see that the two exponents rf^ and ds characterize the important quantities 
(v(t) ■ v(0)) and P{t) for the diffusion process in a fractal medium. 

As p is increased above Pc, the percolation cluster becomes nonfractal, and the particle 
on it starts to diffuse normally at large length scales. In this case, we have a random walk 
with fixed scatterers placed with weakly correlated disorder, and thus we may describe this 
regime as the Lorentz gas regime |TO| , p!T| . In a Lorentz gas, (v(t) • v(0)) generally exhibits a 
long-time tail, obeying the following power law: 

(v(t)-v(0))~r'^/2-\ (4) 

and P{t) scales as in normal diffusion, 

P{t) ~ (5) 

where d is the Euclidean dimension of the cluster. 

Previously diffusion in the two regimes (the anomalous diffusion regime and the Lorentz 
gas regime) were often studied by using very different techniques. We show in this paper, 
the method of spectral analysis of the hopping transition probability matrix provides an 
elegant way of studying both regimes using the same method. 

This work is organized as follows. In Section we review and summarize important 
properties of the matrix W, and in Section ^ we discuss the significance of the eigen- 
spectrum of W. The particular significance of the subdominant eigenvalue Xmax is shown 

4 



in Section |^ In Section 0, the Arnoldi-Saad approximate diagonalization scheme is dis- 
cussed, and the main numerical resuhs of this work are given in Section [V^. A summary 
of our findings is presented in Section |V11| . A prehminary account of this work was given 
where some of the present results were already reported. For completeness, we also 
reproduce those results here with more details, together with new ones particularly on the 
scaling of the subdominant eigenvalue. 



II. HOPPING TRANSITION PROBABILITY MATRIX 

Since the transition probability matrix W contains all the information about the ge- 
ometry of the underlying substrate on which the random walk takes place and also the 
kinematics of the walk itself, all the characteristics of the walk can be extracted from this 
matrix. Because of the importance of W we enumerate in this section some properties of 
W, especially of its eigenvalue spectrum. Much of the information is hidden in its eigenvalue 
spectrum as will be elucidated later. Some of the properties discussed here are already well 
known but are reproduced here for completeness, while others are not commonly appreciated 
despite their potential usefulness. 



A. Basic properties 

By definition, Markov matrices must have the following two properties. First, for all i 
and j, 

Wij > (6) 

and second, 

i 

The second property ensures the presence of the eigenvalue 1 with at least one normalized 
left eigenvector having all components positive and of the same magnitude. There is also 
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a corresponding normalized right eigenvector with eigenvalue 1. This right eigenvector in 
fact represents the equilibrium probability distribution of the diffusing particle. Thus, for 
a connected cluster (with no isolated, unreacheable part), the corresponding eigenvector of 
W can be made positive-definite. 

It is well known that the eigenvalues of the Markov matrix satisfy 

|A| < 1. (8) 

For completeness sake a proof is given here. Let uj be the left eigenvector of W corresponding 
to the eigenvalue A, then 

J2WijUj = \ui. (9) 

j 

Taking the modulus of the above equation, we get, 

l^WijUjl = \\\\ui\, (10) 
j 

J2W,,\u,\ > \X\\u,\. (11) 

j 

Now summing over the index i in the above equation, we get 

EK-lEw^^.>|AlEKI> (12) 

j i i 

3 « 

We have made use of the second property of the Markov matrix. The result Eq. (^ then 
follows immediately. 

B. Detailed and extended detailed balance 

The transition matrix W is said to satisfy the detailed balance condition if its elements 
satisfy the following condition for all (i, j): 
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W,^p^^ = W^,pl (14) 

where p| is the j-th element of the equihbrium probabihty distribution, i.e., the normahzed, 
positive-definite right eigenvector of W for the eigenvalue 1. We now show that W satisfying 
the detailed balance condition can be transformed into a symmetric matrix by a similarity 
transformation, thus proving that all eigenvalues of W are real. 

A proof goes as follows: Taking the square root of the above equation we get 

iw.^Y'\p^)"' = {w,.y/'{pty^'. (15) 

Now, let the matrix A be diagonal with the i-th diagonal element equal to {PiY^'^- Then its 
inverse A~^ is also diagonal with its j-th diagonal element equal to (p|)~^/^. We will show 
that the matrix A~-'^WA is symmetric. Let us consider the ij-ih element, {pf)~^/'^Wij{p'^Y/'^ , 
of this matrix: 

= (f^)-'^'wj,(p^y/^ (16) 

which is just the ji-th element, thus completing the symmetry proof. 

The matrix W is said to satisfy the extented detailed balance condition of order n if it 
satisfies the following property for all (i, j)- 

{wnijp'j = {w^jiPt (17) 

for p = n but not for p = n — 1. Prom the previous argument on the detailed balance 
condition we can say that all eigenvalues of W" are real if W satifies this condition. But 
W will in general have complex eigenvalues in conjugate pairs (as W itself is real). Thus all 
complex eigenvalues of W are such that A" are real, as the eigenvalues of W"' are just the 
n-th power of those of W. Therefore, all the complex eigenvalues of W which have the same 
magnitude converge to either a single real eigenvalue for W" (which is A") or to the two 
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values ±|A"| . Since the eigenvectors of W and W" are the same, all the linearly independent 
eigenvectors corresponding to the different eigenvalues having the same magnitude will have 
the same one or two eigenvalues, producing a degeneracy of order between n/2 and n for 
W". 

The detailed balance condition is satisfied by most models of diffusion (such as the blind 
and myopic ant random walks considered here) , but it is possible to violate it by considering, 
e.g., directed wsXks. A particularly simple example is a random walk taking place on a unit 
square but with a preference for moving clockwise. In this case, the detailed balance is 
not obeyed, but the extended detailed balance is obeyed at the 4th order. More involved 
examples may be found when the kinetic rules of diffusion induce the formation of vortices, 
reminiscent of what actually happens in turbulence. 

C. Oscillations in autocorrelations 

The position autocorrelation function can be expanded in terms of the real (positive and 
negative) and complex eigenvalues in the following way: 

(r(n) • r(0)) = ^ a(A+) A!^ + (-1)" ^ «(A-) A!! 

+ 2a(Ac, uj) cos[ujn + 0(Ac, a;)] A", (18) 

where A+ and — A_ are real, positive and negative eigenvalues respectively and Ace*'^ denotes 
a complex eigenvalue with modulus Ac- The last sum is over the complex eigenvalues and 
thus relevant only for the matrices W which do not satisfy the detailed balance condition. 
The cosine term arises in this case because the complex roots come in conjugate pairs. This 
term, if present, gives rise to oscillations in the position autocorrelation function, different 
from the even-odd oscillation that stems from the negative real eigenvalues. This oscillation 
is persistent (shows up even in the long time limit) if the prefactor of the cosine term A" is 
equal to one. For A < 1, the oscillation is exponentially damped. 

Thus W (if it does not satisfy the detailed balance condition) can have both persistent 
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and transient oscillations in the position autocorrelation function of period larger than two. 
Now, if W satisfies the extended detailed balance of order n, then the autocorrelations will 
have oscillations of period at most 2n, again, either persistent or transient depending on the 
modulus of the eigenvalue causing the oscillation. In this way, the order of the extended 
detailed balance and the period of oscillation in the autocorrelation are closely related. This 
is intuitively acceptable when we consider the case of vortices for a typical example of W 
satisfying the extended detailed balance condition. 

D. Myopic and blind ant 

In light of the above discussion, let us now concentrate on the properties of W for the 
myopic and blind ants specifically. Blind ant in a disordered medium has a symmetric W 
as all the nonzero off-diagonal elements are equal to the same value, 1/z (where z is the 
coordination number of the lattice). Thus all its eigenvalues are real. It does not have 
an eigenvalue equal to —1, so its position autocorrelation function can only have transient, 
even-odd oscillation. 

On the other hand, W for the myopic ant is not symmetric as the hopping probability 
depends on the location of the random walker, that is, it is equal to 1/zi, where Zi is the 
number of occupied nearest neighbors of the particular site. However, it does satisfy the 
detailed balance condition, and thus can be transformed into a symmetric matrix by a 
similarity transformation as seen before. Thus the myopic ant also can have at most an 
even-odd oscillation. However, such an oscillation can be persistent since an eigenvalue — 1 
is possible. This occurs when the cluster is bipartite. This is because the myopic ant on a 
bipartite cluster have eigenvalues in pairs of A and —A, as can be seen easily by converting 
an eigenvector corresponding to eigenvalue A into one corresponding to eigenvalue —A by 
simply making the components corresponding to the sites of one of the subclusters opposite 
in sign. Myopic ants on a nonbipartite cluster, however, can only have a transient, even-odd 
oscillation. 
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III. EIGENSPECTRUM OF W AND DYNAMICS OF THE WALK 



Now we embark on the task of understanding the eigenspectrum of W and building up 
the connection between the spectrum and the characteristic functions of the walk in the time 
domain (such as the velocity autocorrelation function and the return to the starting-point 
probability). 

A. Density of eigenvalues n(A) 

In this work, let us define the density of eigenvalues ^(A) as the number of eigenvalues 
of W normalized by the range of A as well as by the number of sites S in the cluster. 
(This differs from the convention used in ||T^ by a factor of S.) The densities n{X) for the 
blind and myopic ants show some interesting common characteristics. First, at the critical 
dilution p = pc, the eigenvalues accumulate as A — 1, which results in a power law increase 
of n(A) as In A ^ 0. As p is increased from its critical value, the piling up of the eigenvalues 
diminishes and the value of n{X) averaged over small regions of A either flattens (in d = 2) 
or actually decreases (in d = 3) as In A ^ 0. This is illustrated in Fig. |I] where n{X) for the 
myopic ant on a small square lattice is plotted for p = 0.593 ~ Pc and p = 0.8 > Pc- Second, 
the eigenvalue spectrum shows sharp peaks at certain characteristic values of A in a region 
far from the asymptotic limit of A — > 1. This is better illustrated in Fig. |^ for the blind ant 
on a small lattice. These observations are disucssed in more details below. 

First, we recall the connection between the density ^(A) and the probability P{t) that 
the diffusing particle returns to the starting point at time t. The mean probability P{t) can 
be expressed in terms of the corresponding probability Pi(t) for a given starting point i as 

P(t) = ^P,(t)p,(0, (19) 

i 

where po{i) is the probability for the particle to be at site i initially. When the initial 
distribution is the equilibrium one (as we normally assume), this probability is uniformly 
equal to 1/5* (where S is the cluster size) for the blind ant while, for the myopic ant, it 
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fluctuates from site to site according to tlie local coordination number. In either case, Pi{t) 
is given by 

P,(t) = efW*e„ (20) 

where ej is the column vector whose components are all zeros except the i-th component 
which is one and eJ is the corresponding row vector of Cj. 

Below we assume the blind ant for simplicity, but the final result is essentially the same 
for the myopic ant as will be noted. Thus, 

P{t) = TtW'/S = ^'i/S. (21) 

i 

Since —1 is not an eigenvalue for the blind ant and there is a buildup of eigenvalues only 
near A = 1 (for p = pc), and since we are only interested in the long time limit where the 
dominant contribution comes from the eigenvalues near one, we can neglect the contribution 
of the negative eigenvalues. Thus summing the eigenvalues in the above equation only over 
A > and by taking the inverse Laplace transform, we get 

2m Jr ' ^o27r7-oo ' 

= Y.KE-\\^\\)/S, (22) 

A>0 

where Eq. (^1]) has been continued analytically to the complex plane and F is along the 
imaginary axis. 

The quantity on the right is the density of eigenvalues as normalized in E (i.e., in the 
space of In A), but for small E (or for eigenvalues near one), it reduces to n(A). Thus we 
see that P{t) is the Laplace transform of the (normalized) n(A), at least for large t where 
the neglect of A < is unimportant. Although the above discussion has been for the blind 
ant, the calculation for the myopic ant is also very similar. In the latter case, P{t) will be 
bounded by the Laplace transform of n{\) multiplied by two different coefficients. 

Since P{t) for the ants in general decays for large t following a power law 0, 

P{t) ~ t-'^^/^ (23) 
11 



the Laplace relationship implies for both blind and myopic ants, 

n{X) r^\\nX\'^'/^-\ (24) 
Similarly for p > Pc, P{t) decays as 

p{t) ~ (25) 

where d is the Euclidean dimension of the lattice, leading to 

n{X)^\\nXf^-\ (26) 

B. Autocorrelation function 

The position autocorrelation function of a diffusing particle can be expressed in the 
following way: 

(r(t) • r(0)) = X^W*Xo + Y'^W*Yo + Z^WZ^, (27) 

where X"^, Y"^, and Z^ are the row vectors whose components are the x, y, and z components 
of the corresponding sites, respectively, while Xq, Yq, and Zq are the column vectors whose i- 
th components are equal to PiX^, pjYi, and PiZ^, respectively, where p is the initial probability 
distribution for the random walker. 

The vectors X, Y, Z can be expanded in terms of the eigenvectors of the transition 
matrix W. For example, we have 

X^E^aUa, (28) 

A 

where is the expansion coefficient and is the normalized right eigenvector with the 
eigenvalue A. Similarly the vectors Xq, Yq, Zq can be expressed in terms of the normalized 
right eigenvectors of W. For example, we can write 

Xo = EcaUa, (29) 

A 
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where cf is the expansion coefficient. 
Thus Eq. ( p7D can be expressed as 

m ■ r(0)) = Y.iKcl + blcl + hlc\)\\ (30) 

A 

noting that the normahzed eigenvectors are orthogonal to each other. Further denoting 

ax = bx- ca, (31) 
where the vector sign and dot product refer to the coordinate space, we can write 

{r{t)-rm = Y.axX'. (32) 

A 

Though we have calculated the position autocorrelation function here, other correlation 
and autocorrelation functions can also be calculated in a similar way. This includes the mean 
square displacement ((r(t) — r(0))^) and position-velocity cross correlation (r(t) ■ v(t)), as 
well as the velocity autocorrelation function. The latter has been shown ||13[ to be related 
to the eigenspectrum of W in the following way: 

(v(t).v(0)) = -5:a,(A-l)2A*. (33) 

A 

Since the eigenspectrum of the myopic ant on a bipartite cluster contains both the 
eigenvalues 1 and —1, after all the transients have died out, its position autocorrelation 
function behaves as 

(r(^)■r(0))~al + a_l(-l)^ (34) 

That is, the position autocorrelation has an undamped oscillation about the value Oi. If 
the cluster is nonbipartite or if the ant is a blind ant, then there is no eigenvalue —1, and 
thus there will be no undamped oscillations. It can easily be seen from Eq. (|33| ) that the 
corresponding velocity autocorrelation for the myopic ant on a bipartite cluster also has an 
undamped oscillation. However, because of the factor (A — 1)^, this persistent oscillation is 
symmetric about zero as shown in |T^. Again, if the cluster is not bipartite or if the ant is 



a blind ant, then such an undamped oscillation will not exist. 
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Since the long time behavior of autocorrelations for the blind ant and the myopic ant 
on a nonbipartite cluster are (normally) dominated by the eigenspectrum near one, the 
summation over the eigenvalues in Eq. ( |33D can be restricted to only the positive ones. 
For the myopic ant on a bipartite cluster, however, there is an accumulation of eigenvalues 
symmetrically about A = ±1. Thus in this case, the factor (A — 1)^ actually makes the 
eigenvalues near —1 dominate over those near one, and the summation can be replaced by 
the one over only the negative eigenvalues. 

Therefore, in the long time limit, the velocity autocorrelation can be written as 

(v(t)-v(0)) = -5:a,(A-l)2A*, (35) 

AeA 

where A is the set of positive eigenvalues near one for the blind ant and the myopic ant on 
a nonbipartite cluster (Case A) and the set of negative eigenvalues near —1 for the myopic 
ant on a bipartite cluster (Case B). Converting the summation in the above equation into 
an integral, 

(v(t) ■ v(0)) ^ - / ciA7r(A)A*, (36) 

where 7r(A) is 

Tt{\)=n{X)axS{\-lf, (37) 

where n(A) is the density of eigenvalues per cluster site. 

This equation can be further cast in the form of the Laplace transform in the long time 
limit, 

|(v(t)-v(0))|^ / rf«e-"V(A), (38) 

where u = | In A | for Case A (blind ant in general and myopic ant on nonbipartite clus- 
ter), and u = I ln|A|| for Case B (myopic ant on a bipartite cluster). This establishes the 
Laplace relationship between 7r(A) and the velocity autocorrelation. Therefore the power 
law behavior of the velocity autocorrelation 
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|(v(t).v(0))|~r^ (39) 

implies a corresponding power law behavior in vr(A), 

7c{X)r^uy-\ (40) 

with u identified appropriately as above. Thus, by measuring the exponent y — 1 for the 
function vr(A), we can deduce the corresponding asymptotic behavior of the velocity auto- 
correlation function. 

The exponent y appearing in the velocity autocorrelation depends on the two cases 



(Case A, at Pc) , , 

(41) 

ds/2, (Case B, at p^) 
for p = Pc, and 

d/2 + 1, (Case A, p > p^) ^^^^ 
d/2, (Case B, p > pc) 
for p > Pc- In Case A, this exponent characterizes the decay of the autocorrelation which 
is always negative for long times (cage effect), while in Case B, it describes the decrease in 
the envelope amplitude of the even-odd oscillation of the autocorrelation function. 

In the latter case, however, the mean of the autocorrelation function from the two con- 
secutive even-odd steps (rather than the envelope) does have the same decay behavior as in 
the first case. Thus, the behavior of the center line for the oscillating velocity autocorrela- 
tion for the myopic ant on bipartite clusters can be obtained from the eigenspectrum near 
1, allowing it to be included with Case A. In fact, some of the numerical results we present 
later demonstrate that such behavior leads to exactly the same exponent y as in Case A. 



C. Local structure in the eigenvalue spectrum 



The eigenspectrum shows sharp peaks at certain characteristic values of A. The A at 
which these peaks occur depends upon the type of lattice and the kind of ant, indicating 
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that the peaks are structures corresponding to small length scales. For example, the blind 
ant on the square lattice shows strong characteristic peaks at A = 0.5 and 0.75. 

The presence of high degeneracies for these eigenvalues can be explained by looking at 
their corresponding eigenvectors. For example, the eigenvectors for the eigenvalue 0.5 reveal 
that the components corresponding to most of the sites are zero and those with nonzero 
values are arranged in pairs of equal magnitudes but opposite signs. In fact, they involve 
the diagonal arrangement of two sites of a unit square having amplitudes of opposite signs 
with the remaining two sites of the square having zero amplitudes. In some cases, there are 
additional pairs of sites with amplitudes of opposite signs hanging off of such arrangements. 
Similarly, when one inspects the components of the eigenvectors for the eigenvalue 0.75, they 
are seen to involve one, two or three dangling ends which are connected to the main part of 
the cluster through the same site which has a zero amplitude. Again, in some cases, there 
are additional pairs of sites with nonzero amplitudes attached to such arrangements. Some 
examples of these arrangements are illustrated in Fig. |^. 

These local structures connecting a small number of sites are simple enough to be present 
in large numbers for essentially all clusters of appreciable size. Thus they give the sharp 
peaks seen in the eigenspectra as in Fig. |I|, but do not contribute much to the long time 
behaviors of diffusion on the structure. 

IV. SCALING FOR THE SUBDOMINANT EIGENVALUE 

As we have discussed, ai p = pc the eigenvalues of W accumulate near one though the 
eigenvalue 1 itself (if it exists) remains isolated from this build up. This behavior near 
the eigenvalue 1 is mirrored in the region near —1 for the myopic ant on bipartite cluster. 
Thus for any finite cluster there is a sharp drop from the maximum eigenvalue 1 to the 
sub dominant eigenvalue Xmax- This gap decreases to zero as the cluster size increases to 
infinity. 

A simple scahng argument for the cluster size dependence of Xmax can be given as follows. 
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In the Laplace relationship between the time domain quantities and the corresponding func- 
tions of A, the variables t and | lnA| are the conjugate variables. Thus, | lnA|"^ corresponds 
to a time scale. In particular, | In Xmax\^^ corresponds to the time scale associated with the 
eigenmode with the eigenvalue Xmax- In terms of the mapping between the random walk 
and sca/ar elasticity problems [T^, the vibrational energy E' of a mode corresponds to | In A|; 



thus again | lnA|~^ corresponds to a time scale. 

Since the subdominant mode for {R{t)^) (which is actually the dominant one for the 
velocity and acceleration autocorrelations) should reflect the geometry and size of the whole 
cluster, the time scale In Xmax should correspond to the time it takes for a diffusing particle 
to just sample the entire cluster of size S. The latter time scales as S'^^'^^f . Therefore, we 
expect 

|lnA^,,| (43) 



where the exponent d^/dj can be replaced by 2/ds for the critical percolation cluster |]T5 
This way of calculating dg is very efficient since it only requires the accurate determination 
of one eigenvalue (the second largest one) in the entire eigenspectrum. The numerical 
confirmation of this power law is described in the section on the numerical results. 

We now propose a more general scaling relation which includes the region near but 
greater than as follows: 

|lnA„,,| = S-^-/''ffiip-p,)Sn, (44) 

where f{x) is a universal scaling function and the exponent a = 1/ (dfu) is the one appearing 
in the usual percolation cluster size scaling 0. 

In order to satisfy the known asympototic behavior for p > pc, we assume that f{x) 
satisfies the following limiting behavior: 

constant, x — > 
fix) = { (45) 

X — s> oo 
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We see that in the x — hmit we regain the previous relation, and the exponent z can 
be determined by demanding that the x ^ oo hmit leads to the normal diffusive behavior 
where = 2 and df = d, the Euclidean dimension. Thus, in this limit we expect 

|lnA^,,| -S-^/'i (46) 

For this to be true, we must then have 

z = {d^/df -2/d)/a. (47) 

Since the cluster is ordered for p > pc, the scalar elasticity analog would be the phonon 
problem. In this case the energy of the phonons E' ~ A;^ ~ where / is a length-scale of 
the phonon, and since the cluster has the dimensionality d, l^ax ~ S^^'^. Thus, we again see 
that Emin ~ I InAmaj:! ~ S*^^/"^ as above. Numerical results supporting the above finiste size 
scaling relation for | In Xmax \ are also given in the section on numerical results. 



V. ARNOLDI-SAAD APPROXIMATE DIAGONALIZATION METHOD 



In this work we have used the Arnoldi-Saad algorithm in order to approximately 
diagonalize the transition probability matrix W. The main advantage of using this algorithm 
in diagonalizing a matrix is that it allows the very accurate determination of a subset of the 
eigenvalues near the maximum in the spectrum. This is convenient since we are interested 
in the asymptotic long-time behavior of the diffusion process which are determined by the 
eigenvalues near the maximum. In particular, we are mainly interested in the asymptotic 
power law behavior of the density of eigenvalues ^(A) and the function 7r(A) as well as that 
of the subdominant eigenvalue Xmax as discussed above. Both because the thermodynamic 
limit requires a large system and because the eigenvalues approach 1 more closely for larger 
size of the cluster, the accurate determination of the asymptotic behavior demands a large 
cluster, and we need a large number of independent large clusters to obtain true quenched 
disorder averages. While the diagonalization of a large number of large matrices is usually 
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very time consuming and inefficient, the Arnoldi-Saad algorithm affords a valuable reduction 
in the size of the matrices to diagonalize without losing the required information. 

It is interesting to note that the region of the spectrum most accurately obtained by our 
approach is exactly the region of most interest. In previous studies of the eigenspectrum of 
the fractal diffusion problem [IT^, the matrix diagonalized was not W but rather a matrix 



constructed from the equation of motion for the analogous scalar elasticity problem [[15 



For the latter, the important region of the eigenspectrum is near zero (recall the relation 
I In A I ~ cu^ where uo is the angular frequency of vibration for the scalar elasticity problem). 
This is where an accurate determination of the spectrum is most difficult by most methods. 

A large amount of CPU time and memory is saved by using the Arnoldi-Saad method 
which allows the reduction of the dimension of the original matrix without compromising 
the needed accuracy of the results. For instance a 5000x5000 W matrix was reduced by 
this method typically to one of size 300x300, and the eigenvalues and eigenvectors obtained 
were reliable to better than one part in 10^. This operation required about 560 seconds of 
CPU time per cluster on a Kubota Pacific Titan P3 mini-super computer or less than 30 
seconds on a CPU of a Cray 2 super computer, including the cluster generation itself. 

The algorithm was described in the literature |T^, but we will sketch it here for com- 
pleteness. In this algorithm, one starts with an arbitrary normalized vector Ui of dimension 
S where W has the dimension of 5* x S*. Then one chooses the subspace dimension m < S 
which is to be the dimension of the reduced upper Hessenberg matrix H which will then be 
diagonalized essentially exactly. The reduced matrix H is obtained recursively and at the 
same time as a sequence of normalized basis vectors Uj (i = 2, m) in the following way: 

First, the element Hu is obtained as ufWui. Second, the product V2 = -ff2iU2 is 
calculated as V2 = Wui — HnUi. We then choose H21 = ||v2|| to satisfy the normalization 
requirement ||u2||^ = 1. At the next iteration, we ffist obtain H12 and H22 as ufWu2 and 
u^Wu2, respectively, and then, U3 and H32 are obtained similarly to the ffist iteration. This 
process continues until all elements of H and all of the Uj are obtained. In general, we have 
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= ufWu,, (48) 

i 

ifj+ijuj+i = Wuj - HijUi, (49) 



where Uj {i = 1, 2, m) form an orthonormal set. 



This procedure ensures that the eigenvalues of H are approximate eigenvalues of the 



original matrix W and that if y is the eigenvector of H with the eigenvalue A, then the 
vector z defined by its components Zj = YnLiluiljUi satisfies the orthogonality 

(W - AI)z • Ui = 0, i = l,2,...,m. (50) 

This means that z is an approximate eigenvector of W to within the subspace spanned by 
the Uj's with the approximate eigenvalue A. 

The approximation is better the larger m is, and also the outer part of the spectrum 



yields better accuracy [|l6l. In practice, the eigenvalues pile up at very high density near the 
outer edge of the spectrum (for pc), which requires with a large (typically greater 
than 100) to be used in place of W to separate the eigenvalues. The accuracy of the 
eigenvalues and eigenvectors can be checked by a number of criteria including one which 
uses the comparisons between the results obtained by using different dimensions m for the 
submatrix H, and also by using different A^. 



VI. NUMERICAL RESULTS 

In this section, we summarize the numerical results for the density of states per site n(A), 
the function 7r(A), and the scaling of the subdominant eigenvalue InAmax- 



A. Density of eigenvalues n(A) 

The density of states n{X) for the blind and myopic ants has been calculated in two di- 
mensions on the square lattice and in three dimensions on the simple cubic lattice, both with 
and without periodic boundary conditions. In practice, n{X) is computed by accumulating 
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the eigenvalues in logarithmic bins over a large number of independent cluster realizations 
and then dividing by the number of clusters, the bin width, and the cluster size. 

In Fig. ^, n{X) at pc is plotted against | lnA| for the blind ant on nonperiodic clusters in 
d = 2 and 3. The averages over 5000, 250 and 900 clusters were taken for the cluster sizes 
S = 400, 800 and 5000, respectively in = 2, while for d = 3 the averages were taken over 
250, 300 and 449 clusters for the cluster sizes S = 400, 800 and 5000, respectively. It is clear 
from the figure that the data for different cluster sizes collapse into a single curve for both 
d = 2 and 3, showing that finite size effects for n{X) are not important. The approximate pc 
used here is p = 0.593 for the square lattice in d = 2 and 0.312 for the simple cubic lattice 
in d = 3. 

In order to evaluate the exponent dg in d = 2, a log-log plot is produced for n(A) as 
a function of | In A| in Fig. ^ for the blind and myopic ants on the square lattice. As seen 
clearly, the data for the blind ant on the periodic and nonperiodic clusters collapse within the 
statistical errors on the same line of slope — 0.35±0.01, showing that boundary conditions do 
not affect the exponent dg. The data for the periodic cluster have been collected on 100 x 100 
grids and the average was taken over 400 clusters and then normalized by the average size 
of the cluster. For the nonperiodic cluster, the average was taken over 900 clusters of 5000 
sites. The data for the myopic ant on the periodic cluster, on the other hand, yield a least 
squares fitted line of a slightly larger slope —0.37 ± 0.01. The dashed lines showing these 
fits had the correlation coefficients greater than 0.999. The data for the myopic ant were 
obtained on 100x100 grids and the average was taken over 250 clusters. We believe that 
the small difference in the slopes between the blind and myopic ants is a reflection of the 
inexact Laplace relationship between n(A) and P{t) for the case of the myopic ant. The 
details of the exponent estimates extracted from these data along with the values of the 
exponent quoted in a standard reference ^ are tabulated in Table |. The blind and myopic 
ants on triangular lattice also had a similar slope but they are not shown in the figure to 
avoid overcrowding. 

The exponent dg for d = 3 has similarly been extracted from log-log plots of n(A) versus 
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I In A I for the blind ant on the periodic and nonperiodic simple cubic lattices and the myopic 
ant on the periodic simple cubic lattice, as shown in Fig. ^. The data for the blind ant on 
the periodic lattice were collected on 30^ grids and the average was taken over 449 clusters. 
The data for the nonperiodic blind ant were collected for 5000 site clusters and the average 
was taken over 500 clusters. The data for the myopic ant were for a 30'^ lattice, the average 
being taken over 500 clusters. The dashed lines shown in the figure are the least square 
fits of the data points and all the fits had correlation a coefficient of greater than 0.999. 
The least square fitted lines for the blind ant on the periodic and nonperiodic lattices have 
the same slope of —0.35 ± 0.01 within the statistical error, showing that the boundary 
conditions do not affect the value of dg itl d = 3 either. Here, however, n(A) normalized by 
the average cluster size for the periodic case does not seem to have the same amplitude as 
the nonperiodic case. The myopic ant on the simple cubic periodic lattice gave a slightly 
larger slope of —0.37 ± 0.01, similarly to the case of d = 2. The details of the extracted 
exponent estimates along with the values of the exponent obtained from a standard reference 
1^ are given in Table |. When these results are compared, we see that they are consistent 
with each other within the respective errors. 

In Fig. 1^, a log-log plot of ^(A) against In A for p > is shown for both d = 2 and 
d = 3. In this case, the substrate becomes more ordered and the n(A) is expected to scale 
as I InAI'^/^"^, where d is the Euclidean dimension of the substrate. Our data show a spiky 
structure {not statistical fluctuations) due to the discrete nature of the eigenvalues. In the 
d = 2 case, the data were collected for p = 0.75 and for p = 0.9 for the clusters on a 100 x 100 
periodic square grid and the average was taken over 250 and 200 clusters. In d = 3, the 
scaling of ^(A) was tested for p = 0.5 and the average was taken over 100 clusters on a 30^ 
periodic grid. In these figures, lines of expected slope (zero for d = 2 and 1/2 for d = 3) are 
drawn to guide the eye. As is evident from Fig. |^, the data for d = 3 have a larger slope 
than for d = 2 as expected. 
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B. Asymptotic behavior of 7r(A) 



As discussed earlier, the function 7r(A) at pc is expected to follow the power law 7r(A) ~ 
I In A|^~^/'^'", where is the walk dimension. This exponent is thus extracted from the slope 
of the log-log plot of vr(A) versus | lnA| as shown in Fig. ^. Plotted in this figure are the 
results from the blind ant for nonperiodic clusters on the square and simple cubic lattices 
and those from the myopic ant for periodic clusters on the same two lattices (for A near 1 
and not -1, as discussed in Sec. [m| ). The dashed lines in the figure are the least squares 
fit and the computed correlation coefficient was greater than 0.999 in every case except for 
the myopic ant in three dimensions where it was 0.997. Some data points toward both ends 
of |lnA| were not included where they correspond to partially filled bins (i.e., only some 
cluster realizations have eigenvalues falling in these bins). The nonperiodic clusters were 
of size 5000 and, in the periodic case, grids of size 100x100 in c? = 2 and of 30^ in ci = 3 
were used. The average was taken over the same number of clusters as in the corresponding 
case of ^(A) except for the myopic ant for the periodic simple cubic lattice where only 200 
clusters were averaged. 

The estimated value oi 1 — 2/dw for d = 2 is 0.30 ±0.01 for the myopic ant on the periodic 
square grid and 0.27 ±0.01 in the case of the blind ant on the nonperiodic square grid. (The 
estimates for the periodic lattice are actually identical for both the blind and myopic ants.) 
In d = 3 the slope of the least squares fit was found to be 0.46 ± 0.02 for the myopic ant on 
the nonperiodic simple cubic grid and 0.43 ±0.01 for the blind ant for the nonperiodic case. 
The exponents calculated from the slopes for the data for the periodic lattice agree very well 
with the previously available ones from Ref. within the respective errors. Although our 
slope estimates themselves are substantially different between the periodic and non-periodic 
clusters, the corresponding exponent d^ calculated from these slopes do not differ quite as 
much. While the remaining deviations are still larger than the statistical errors and thus a 
source of concern, we believe that the periodic lattice results are inherently more trustworthy 
and the non-periodic lattice results are somehow refiecting the more severe boundary effects 
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and/or the effects of the cluster geometry imposed by arbitrarily stopping the cluster growing 
algorithm at a predetermined size. We do not, however, understand why this effect seems 
to be much more prominent for vr(A) than for n{X). 

In Fig. I a log-log plot of 7r(A) against | lnA| for p > pc is presented. For p > pc, the 
function 7r(A) is expected to behave like | InAj'^/^ as A — 1. Again in this case we find some 
structure in the data, and lines of slope 1 and 1.5 are drawn in the figure to guide the eye. 
For d = 2, the data were collected at p = 0.75 and p = 0.9, and for d = 3, the data were 
collected for p = 0.5. The number of clusters over which the average was taken is the same 
as the corresponding data for n{X). Again it is evident from the figure that the slope of 7r(A) 
in the case of c? = 3 is larger than in c? = 2 as expected. The summary of all these results 
are given in Table |I[ 

C. Subdominant eigenvalue Xmax 

The values of the maximum eigenvalue below one, Xmax, in two and three dimensions are 
plotted in Fig. ^ against the cluster size. They were calculated for clusters of size 100, 200, 
400, 800, 1000, and 5000 for the myopic ant on the triangular lattice and for the blind ant 
on the square lattice for cluster sizes of 100, 200, 400, 800, and 5000. In three dimensions, 
Xmax was found for the blind ant on the simple cubic lattice of cluster sizes of 100, 400, 
1000, and 5000. The averages for the cluster sizes 100, 200, 400, 800, 1000, and 5000 were 
taken over 10000, 5000, 5000, 2500, 2500, and 1000 clusters, respectively. The accuracy in 
the numerical evaluation was better than one part in 10^ in all cases. 

As discussed above the subdominant eigenvalue In A scales with the size of the cluster, 
S, as I InAmaxI ~ S~'^'"^'^f dX p = pc- We have extracted this exponent d^/df from the slope 
of the log-log plot of I In Xmax I versus S as shown in Fig. |1^. The straight lines drawn are 
linear least square fits and all of them have a correlation coefficient of greater than 0.9999. 
The statistical errors are smaller than the symbol sizes in all cases. 

In (i = 2 we have obtained a slope magnitude of 1.546 ± 0.019 for the myopic ant on the 
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triangular lattice and a slope of 1.529 ± 0.004 for the blind ant on the square lattice. In 
comparison, using the value of dg = 1.30±0.002 as given in Ref. and assuming the scaling 
relation ds = 2df/dw, we would get the corresponding value of d^/df = 1.538 ±0.024. Thus 
we see that our estimate computed from the finite size scaling of the subdominant eigenvalues 
is in excellent agreement with the direct simulation result (and the scaling relation used). 

In d = 3 the value of the slope magnitude obtained is 1.508±0.007. Again using the value 
of dg = 1.328 ± 0.006 from Ref. and the relation dg = 2df/dy^, the simulation estimate 
of this quantity would come out to be 1.506 ± 0.007, again in excellent agreement with our 



analysis. The summary of these results are given in Table |ITI . 

In order to check the finite size scaling relation for | In Xmax I for p close to but greater 
than pc, we plot in Fig. [l^, | \n\max\S'^^'^'' against {p — pc)S" for p = 0.64, 0.65, 0.66 in d = 2 
and for p = 0.33, 0.34, 0.35 in d = 3, where a = l/{dfv). All the data in d = 2 are for 
the blind ant on square, nonperiodic clusters, and in d = ?>, for simple cubic, nonperiodic 
clusters. In d = 2, the average was taken over 10000, 10000, 5000 and 1000 clusters for 
cluster sizes of 100, 400, 1000 and 5000, respectively. In = 3, the average was taken over 
100, 250, 300 and 1000 clusters for the cluster sizes of 100, 400, 1000 and 5000, respectively. 
The statistical errors are typically of the order 10~^ and thus less than the symbol size in 
each case. 

As is evident from the figures, both in d = 2 and 3, the data for different values of 
p collapse onto a single curve, which is the universal scaling function f{x) (where x = 
{p—Pc)S'^), supporting the scaling form given by Eq. (0). The scaling function is numerically 
fitted to a quadratic curve of the form ax'^ + bx + c as shown by the dashed lines in the figure. 
The values for the coefficients a, b, c have been numerically estimated as 2.646, 0.893, 2.358 
from the square lattice data and 0.894, 0.276, and 1.044 from those on the simple cubic 
lattice, respectively. Even though f{x) is expected to depend only on the dimensionality, 
in general there will be lattice-dependent metrical factors which will affect the values of the 
coefficients. 

The nonzero value of c indicates that f{x) indeed tends to a nonzero constant in the limit 
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X ^ 0. This is consistent with the scahng relationship for | In Amaa l at j9 = pc- Our scahng 
prediction for the value of the exponent z is 1.21 for d = 2 and 1.81 for d = 3, respectively, 
from Eq. (^Sf). The value of z was not estimated numerically since it is valid only in the 
asymptotic limit a; ^ oo thus needing the calculation of Xmax for cluster sizes larger than 
available due to computational limitations. 



VII. SUMMARY 



In summary, we used the Arnoldi-Saad algorithm [0] to approximately diagonalize the 
hopping transition matrix and successfully extracted the dynamical exponents rf^ and ds at 
the critical value of p for the percolation cluster. We found that, for p > pc, the eigenvalue 
spectrum reflects the discrete, local lattice structure more clearly, and mainly because of 
this there is a spiky nonstatistical structure in the density of states. Despite the numer- 
ical difficulty caused by this, we clearly observe the crossover between the fracton regime 
(anomalous diffusion) and the phonon regime (Euclidean diffusion with a long time tail). 

We have also shown that the finite size scaling of the sub dominant eigenvalue Xmax to 
be the most computationally efficient method to extract the critical exponent combination 
dw/df (which is equal to 2/ds for percolation). The computational efficiency is due to the 
fact that only the largest eigenvalue (below one) has to be evaluated. 

The method of eigenspectrum analysis as presented here reduces the study of diffusion 
in complex, random media to that of a numerical analysis of large, random matrices. Such 
an analysis can be carried out with very high precision in a systematic way taking the 
best advantage of modern computing equipment (with vector and parallel processing). In 
fact, this approach was recently applied by us to a number of problems including the scalar 
elasticity of loop- enhanced stTuctnies and a study of the breakdown of the scaling relation 
ds = 2df/duj for treelike structures |]19[. Further applications contemplated include the 



problem of diffusion in an external velocity field, a possible model of Brownian transport in 
the presence of macroscopic convection. 
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APPENDIX: EXPANSION COEFFICIENTS ax FOR PERIODIC BOUNDARIES 



We show in this Appendix how the calculation of the expansion coefficients a\ which 
appear in the function vr(A) can be obtained from the velocity vectors in the case of the 
periodic boundary condition. From Eqs. (^) - (|3TD , we have 

OA = E vxis)r{s) ■ E ux{so)p{so)r{so), (Al) 

s so 

where ua and va are normalized left and right eigenvectors of W for the eigenvalue A, p 
is the initial probabihty distribution, and the dot product refers to the coordinate space as 
before. 

Replacing the coordinate vectors r(s) by the velocity vectors [r(s) — r(sj)] where the sites 
Si are the nearest neighbors of the site s, the ffist of the two factors in the above equation 
becomes 

s s ^\^) i=l s ^{^J i=l 

Here z{s) is the number of occupied neighbors of the site s in the case of the myopic ant. 
For the blind ant, z{s) stands for the lattice coordination number and, if the ith neighbor 
is not a cluster site, then Sj = s is understood. 

The second term on the right can be simplified by exchanging the summations over s 
and sf. 

E -a(^)-^ E = E E Z^M^'^ 

s ^y^J i=l s' i=l '^V^i) 

= Y.<s)\vx{s') 

s' 

= XY.r{s')vx{s% (A3) 

s' 

where we used the relation 

i = l s" 

Therefore we get 

28 



E^A(^)r(.) = 7^^E^A(.)^E[r(^) (A5) 

s U - Aj g Z[S) .^-^ 

Similarly the second factor in Eq. is 

T.Ms)r{s)p{s) = -^Y.Ms)pis)^^Y.l^^^^ (A6) 

s U - Aj ^ z[s) -^^ 



Once in this form, the absolute coordinates are not required; rather, only the nearest neigh- 
bor velocities are. In the calculation of the expansion coefficients ax, then, we only need to 
keep track of the velocity vectors in the case of the periodic boundary conditions. 
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FIGURES 

FIG. 1. Density of eigenvalues n(A) for the myopic ant on the square lattice cluster of size 200 
at p = 0.593 « pc and p = 0.8. This figure shows the build up of n(A) near A = ±1 at p = 0.593 
and the flattening of n(A) for p > pc- 

FIG. 2. Density of eigenvalues n(A) at p = 0.593 for the blind ant on small square lattice 
clusters of size 800 and 100 sites. This figure shows the peaks due to the local structures of the 
cluster as well as the increase in n(A) near one for larger size clusters. 

FIG. 3. Illustration of the eigenvectors which contribute to the peaks in the eigenvalue spectrum 
at p = 0.75 and p = 0.5 using a 100-site cluster on the square lattice (sites shown with dots). Some 
(not all) such eigenvectors are chosen and the sites of nonzero amplitudes for eigenvectors with 
eigenvalue A = 0.75 are marked by the symbol + and those of A = 0.5 are marked by A. 

FIG. 4. Density of eigenvalues n{\) normalized by the cluster size for the blind ant on nonpe- 
riodic clusters on square and simple cubic lattices. The upper curve corresponds to d = 3 and the 
lower one to d = 2. The symbols +, A and o correspond to the cluster sizes 400, 800 and 5000, 
respectively. 

FIG. 5. Log-log plot of ra(A) against | In A| in d = 2 at 7? = 0.593 ^ Pc- The symbols A 
and + correspond to the blind ant on periodic and nonperiodic square lattices and the symbol o 
corresponds to the myopic ant on periodic square lattice. 

FIG. 6. Log-log plot of n(A) against | lnA| in d = 3 at p = 0.312 pc- The symbols A and 
-|- correspond to the blind ant on periodic and nonperiodic simple cubic lattices and the symbol o 
corresponds to the myopic ant on periodic simple cubic lattice. 

FIG. 7. Log-log plot of n(A) against |lnA| in d = 2 and 3 for p > pc- The symbols o and A 
stand for the blind ant on the periodic square lattice aX p = 0.75 and p = 0.9, respectively, and 
the symbol -|- stands for the blind ant on the periodic simple cubic lattice at p = 0.5. The lines of 
slopes zero and 1/2 are drawn to guide the eye. 
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FIG. 8. Log-log plot of 7r(A) against |lnA| in d = 2 and 3 at ^? = pc- The symbols A and 
o stands for the myopic ant on the periodic square lattice and the blind ant on the nonperiodic 
square lattice, respectively, and the symbols x and + correspond to the myopic ant on the periodic 
simple cubic lattice and the blind ant on the nonperiodic simple cubic lattice, respectively. 

FIG. 9. Log-log plot of 7r(A) against | In A| in d = 2 and 3 for p > pc- The symbols o and A 
stand for the blind ant on the periodic square lattice aX p = 0.75 and p = 0.9, respectively. The 

symbol + stands for the blind ant on the periodic simple cubic lattice at p = 0.5. 

FIG. 10. Log-log plot of I In A 

max\ against the size S of the cluster. The symbol o corresponds 
to the myopic ant on the triangular lattice. The symbols A and -|- correspond to the blind ant on 
the square and simple cubic lattices, respectively. 

FIG. 11. Finite size scaling of |lnAj7iaa;| with respect to the size of the cluster and {p — Pc)- 
Quadratic fits for the universal scaling functions f{x) are shown in dashed lines, where 
X = (j) — PcjS"^ and the upper curve corresponds to d = 2, the lower to d = 3. 
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TABLES 



TABLE I. Estimates of dgjl — 1 from the density of eigenvalues n(A) 



d 


This work 


Using dg from Q 


d/2 - 1 (for p > pc) 


2 


-0.35 ±0.01 


-0.347 ±0.013 





3 


-0.35 ±0.01 


-0.336 ± 0.003 


1/2 




TABLE IL 


Estimates of 1 — 2/d^ from vr(A) 


d 


This work 


Using dw from J 


1 d/2 (for p > Pc) 


2 


0.30 ±0.01 


0.30 ± 0.01 


1 


3 


0.46 ±0.01 


0.46 ± 0.03 


3/2 




TABLE III, 


. Estimates oi d^/dj from \n\max 


d 


Ant type 


This work 


Using 2/ds from Q 


2 


myopic 


1.546 ±0.019 


1.538 ±0.024 


2 


bUnd 


1.529 ±0.004 


1.538 ±0.024 


3 


blind 


1.508 ±0.007 


1.506 ± 0.007 
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